Prom String Nets to Nonabelions 
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We discuss Hilbert spaces spanned by the set of string nets, i.e. trivalent graphs, on a lattice. 
We suggest some routes by which such a Hilbert space could be the low-energy subspace of a model 
of quantum spins on a lattice with short-ranged interactions. We then explain conditions which 
a Hamiltonian acting on this string net Hilbert space must satisfy in order for its ground state 
and low-lying quasiparticle excitations to be in the DFib topological phase. Using the string net 
wavefunction, we describe the properties of this phase. Our discussion is informed by mappings of 
string net wavefunctions to the chromatic polynomial and the Potts model. 



I. INTRODUCTION 

The mathematical theory of anyons - modular tensor 
categories - is extremely rich in examples, and existing 
physical theory provides a way of describing most of these 
examples as effective Chern- Simons gauge theories. In 
contrast, our knowledge is limited when it comes to iden- 
tifying plausible solid state Hamiltonians from which a 
(2D) state of matter could emerge whose effective low 
energy description is, in fact, a Chern-Simons theory. 
The off-diagonal conductivity of fractional quantum Hall 
(FQHE) systems is tantamount to the equation of mo- 
tion for a Chern-Simons Lagrangian, so Hall systems are 
the most developed source of such examples. The best- 
studied example among abelian states is the v = 1/3 
Laughlin state. The foremost candidate among non- 
Abelian states is the Pfaffian state 3 -, which is believed 4 ^ 
to be realized at the v = 5/2 fractional quantum hall 
plateau. Beyond v — 5/2, more delicate plateaus at 
v = 12/5,4/7, etc. may also support nonabelions. How- 
ever, in this paper we explore a quite distinct family of 
Hamiltonians. 

Because magnetic interactions in solids can be at en- 
ergy scales as high as ~ 10 3 Kelvin, it would be very ex- 
citing to find realistic families of spin Hamiltonians rep- 
resenting a nonabelian phase. (This has essentially been 
accomplished^ for the simplest abelian phase, gauge 
theory, although the corresponding experimental system 
has not been clearly identified.) This goal has been pur- 
sued for several years through the study of model Hamil- 
tonians H acting in an effective Hilbert space H whose 
degrees of freedom are either unoriented loopsi, or, more 
recently, branching loops called "string nets'—. Such 
Hilbert spaces Ti. are a kind of half-way house. Even- 
tually, it will be necessary to understand how local spins 
can encode effective loops and nets, and some ideas on 
encoding nets are presented in Section II. However, the 
premise of this paper is that we already have a Hilbert 
space Ti. spanned by the simplest type of string net G - 
the lines are unoriented and unlabeled, the nodes have 
valence 3 and lack internal states. Our goal then is to 
formulate, in the most general terms, what properties a 



Hamiltonian H : Ti — ► Ti should have to describe the 
simplest topological phase of string nets, the "doubled 
Fibonacci theory", DFib. DFib is not only nonabelian 
but actually computationally universal^ and thus an ex- 
tremely attractive target phase. The approach here is 
intended to complement the ideas presented in Refs. 
and [J 

To summarize our approach in a phrase: "nature ab- 
hors a degeneracy" . (Consider, for example, eigenvalue 
repulsion for random Hamiltonians). There is an irony 
here because topological phases are nothing else than a 
degenerate, yet stable, ground state for which no classical 
symmetry exists to be broken. From this viewpoint, we 
will see that building DFib (and other phases?) amounts 
to setting a trap for nature. By compelling a certain 
space V{D, n) of low energy modes for a given fixed 
boundary condition (specified by integer n, as we de- 
scribe below) to have dimension equal to d, unexpectedly 
small, we trap a class of Hamiltonians into an exponential 
growth of degeneracy: lim n ^ 0o (dim(V'(.D,n))) 1 ' n = t = 
■ We present a rather surprising derivation of DFib 
from dimensional considerations alone; the 6j-symbol de- 
rives from the assumption of unitarity and minimal di- 
mension of "disk spaces" . Although the 6j-symbol obeys 
the pentagon equations we do not use the pentagon equa- 
tion to find the 6j-symbol. The physical significance is 
that DFib should be a robust phase stabilized by a type 
of eigenvalue repulsion. 

After DFib is derived, a beautiful formula of Tutte 
(compare Ref. 0) allows us to make an exact connec- 
tion to the Q = t + 2 w 3.618 state Potts model, where 
t = i * ne golden ratio. We find that the exactly 
solvable point in the DFib phase is the high temperature 
limit of the low temperature expansion of the r + 2-state 
Potts model. This allows us to conjecture that the topo- 
logical phase extends downward in "temperature" until 
the critical point is reached at log (3 = \Jt + 2 + 1. This 
suggests a one-parameter family of DFib-Hamiltonians 
whose ground state wave functions are not strictly topo- 
logical but have a "length" or bond-fugacity, x, satisfying 
0.345 w 1/(a/t + 2 + 1) < x < 1, implying a considerable 
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stability within this phase. 

The paper is organized as follows. In Section II, some 
ideas are presented for how string nets could emerge from 
microscopic models of quantum spins on a lattice. In Sec- 
tion III, DFib is derived from dimensional considerations. 
In Section IV, we describe, in string net language, the 
quasiparticle excitations of DFib. In Section V, string 
net wavefunctions and their squares are related to the 
chromatic polynomial. In Section VI we show that the 
topological string net wave function has a "plasma anal- 
ogy" to the T+2-state Potts model, and in Section VII we 
discuss our conclusions. In the appendix, Baxter's hard 
hexagon model is used to extend a theorem of Tutte's. 



II. HOW TO CONSTRUCT A NET HILBERT 
SPACE n 

In this paper, we will be concerned with Hamiltonians 
H acting on Hilbert spaces Ti (£) of wave functions 4* that 
assign complex- valued amplitudes to string nets ("nets" 
for short) on a surface S. The surfaces £ which could be 
relevant to experimental systems are presumably planes 
with some number of punctures. However, it is quite 
profitable conceptually and for the purpose of numerical 
simulations to think about higher-genus surfaces as well. 
A net is what mathematicians call a trivalent graph; it 
has only simple branching and no "dead ends" (univalent 
vertices) except as defined by boundary conditions at the 
edge of the surface. According to this definition, (the 
Dirac function on the net in) Figure ^ a) is not in the 
Hilbert space Ti but (those of) Figure Q]b) and He) are. 
In most of this paper, we will simply assume that the 



are sketched for writing a spin Hamiltonian K : H — > 
H on a large Hilbert space H of microscopic degrees of 
freedom so that the ground state manifold Ti of K will 
be the "Hilbert space of nets" on which this paper is 
predicated. 

All our K break S'[/(2)-invariance and require fine tun- 
ing. Ideas A and B are conceptually very simple but both 
require a three-body interaction. Idea C is really an en- 
cryption of B into a 2-body interaction on higher spin 
(spin = 3/2) particles. 



H = <8>bondsC 2 , i.e. is a Hilbert space of spin = 1/2 
particles living on the links of a trivalent graph such as 
the honeycomb. We interpret an Sf = 1/2 link as one 
on which the net lies. We take K = T, s i tes K s where K s 
projects onto the subspace of C 2 ® C 2 ® C 2 (i.e. of the 
Hilbert space of the three spins surrounding a lattice) of 
total J2i=i^i = — 1/2- This forbids dead ends. Unfortu- 
nately, the function 



total S z 



eigenvalue 
-3/2 
-1/2 
1/2 
3/2 



energy 
0, 

nonzero 







(1) 



is not quadratic (no parabola passes through 
(-3/2,0), (-1/2, nonzero), (1/2,0), and (3/2,0). Hence, 
when K s is expanded in products of a z i, i — 1,2,3, 
running over the bonds meeting s, it must contain a 
cubic cr z ia z 2cr z 3 term. 



a) b) c) 




FIG. 1: a) is not in the Hilbert space Ti, but b) and c) are. 
In these figures, the outermost circle is the boundary of the 
system, where nets are allowed to terminate. The endpoint 
in the middle in (a) is a violation of the "no dead ends" con- 
dition. 

Hilbert spaces Ti(T,) arise as the low-energy subspaces of 
the Hilbert spaces of a system of spins or electrons in a 
solid or ultra-cold atoms on an optical lattice. In such 
a formulation Figure [I] a) could be thought of as a high, 
but finite, energy state of the spins, electrons in a solid 
or ultra-cold atoms. 

In this section, however, we will consider the question: 
from what kinds of lattice models do nets emerge in the 
low-energy description so that the Hilbert spaces W(S) 
are the low-energy subspaces? Three ideas A, B, and C 



B. 

An alternative is to put a spin = 1/2 particle at the 
sites of the honeycomb, so Ti = ® s itcsC 2 , C 2 = (+, — ). 
For each consecutive pair p of bonds, K has a term K p : 
K = "EpKp, where K p is a diagonal matrix all of whose 
entries are except for that corresponding to the illegal 
pair of edges shown in fig. 0). K penalizes both isolated 




FIG. 2: An illegal pair of bonds which is energetically penal- 
ized by K p . 

+'s and +'s with exactly one + neighbor. The latter 
situation is shown in fig. [3] Interpreting those bonds 
bounded by two + signs as the ones on which the net lies, 
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FIG. 3: An up-spin which has only a single up-spin neighbor. 
This is energetically penalized by K — Y, P K P . 



we see that the zero modes of K are precisely the nets 
(trivalent graphs) within the honeycomb. Unfortunately 
K seems resolutely 3-body. 



s q t r 

-3/2 -1/2 +1/2 +3/2 
FIG. 5: the bottom row represents the spin eigenvalues 



B'. 

To set the stage for our final construction, it is help- 
ful to reverse + and — spins on the index 2 Bra- 
vais lattice L' within the honeycomb, honeycomb = 
L U V. With this convention, if = £ p contercd on L K p + 

Dp- centered o„ v where projects to the highest 
total S z eigenvalue, S z = 3/2, and K' p > projects to the 
lowest total S z eigenvalue, S z = —3/2. In other words, 
the penalized configurations are 3 consecutive pluses or 3 
consecutive minuses. Also, we now draw bonds between 
plus sites of L and minus sites of L' (L and L' labeled 
as black and red respectively in Figure |2J. K is still nec- 
essarily 3-body but at least it now has the form of the 
Klein Hamiltonian, see e.g. ref. 




FIG. 4: Upon reversing the spins on one sublattice of the 
honeycomb lattice, K p and K' p i now penalize maximum and 
minimum S z eigenvalues respectively. 



c. 

The idea here is to take adjacent pairs of site spins 
from B' and encrypt them as the state of a spin = 3/2 
particle living on the bond b joining the two sites. We 
orient b from L' to L. Let us set up an indeterminate 
bijection (fig. |5j). We now "simulate" K on a Hilbcrt 
space ft = ®bond s C 4 , the space of a spin — 3/2 particle 
on each bond. On sites I € L (/' S L') we must penalize 
t®t (q®q)- 



Furthermore we must penalize inconsistent encryp- 
tions. For I G L (I' S L') the following pairs yield incon- 
sistent site labels: q® r, r ® q, q® t, t ® q, r, r ® s, s®t, 
and t <£> s (q <£> s, s ® q, q (£> t, t (£> q, r ® s, s ® r, r ® i, and 
f <g> r). 

We again obtain the space of nets as zero modes H CH 
by fixing if to be the following 2-body Hamiltonian: 

k = (lit®* + n g ® r + n r( g„j + n^^i + n t ^, 9 

+ n s( g, r + n r( g, s + iLjgi + n t S ) 

+ (n g( gig + n 9 ^, s + n s ® q + n g ^, t + n t( g, g 

I'eL' 

+ n r gi s + iis^j. + n r t + n t r ), (2) 

where IL x tg>y denotes the projector onto the state x® y. 



spin j 




-3/2 -1/2 1/2 3/2 

spin i 



FIG. 6: The horizontal axis labels the ith spin and the vertical 
axis the j'th spin. 

We can write a Hamiltonian which effectively accom- 
plishes such a projection in terms of the spins Sf (i is the 
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bond index). For black sites H has terms 

((s? + i) 2 + (s; + i) 2 -i/2) 

(Sf + S| - 2) (Sf + S ZJ - 3) (3) 
and for red sites 

((sf-i) 2 + (s;-i) 2 -i/2) 

(<?*+<?* + 2) (Sf+S* +3). (4) 
The origin of these terms is illustrated in Figure 

III. H, H, AND V 

From now on, we will be concerned with he proper- 
ties of the Hilbert spaces H(E) on a surface E. If E has 
boundary, <9E, we fix a boundary condition by specifying 
points where the nets must end. In the case where S has 
connected boundary and the boundary condition consists 
of n points, we denote the Hilbert space by 7i(E, n). We 
will be interested in "isotopy invariant" wave functions 
whose value is independent of deformation. (In Sec- 

t((D) = ?(©) 
t(®) = 

FIG. 7: The condition of isotopy invariance. Wavefunctions 
assign a complex amplitude for any string net (which is the 
amplitude for this configuration to occur). These equations 
mean that wavefunctions in the low-energy Hilbert space as- 
sign the same amplitude to two string nets if one can be ob- 
tained from the other by smooth deformations. 

tion VI we relax this condition to allow certain bond fu- 
gacities.) To avoid unnormalizable wave functions, these 
loops should really live on a lattice, as in the previous 
section. 

To produce invariant VP's we consider Hamiltonians 
H : 7i — > 7i which contain fluctuations sufficient to en- 
force isotopy invariance (fig. EJ on all low energy states. 
It should be remarked that it is not easy to set up such 
terms on a lattice; some fine tuning may be required. 
Also there are questions of ergodicity - H must have suf- 
ficient fluctuations that crystals do not compete with the 
liquid condition described by Fig. [7] Nevertheless we 
start by assuming these problems solved: that we have 
7i and a family {H} whose ground states V consist of 
isotopy invariant wave functions 'J. 

The "axioms" we impose on the Hamiltonian H are 
implicit in the following conditions that we require its 
ground state manifolds V(E, n) C H(T,, n) to satisfy. 

Axiom 1. H is gapped - this makes F(E, n) sharply 
defined. 



Axiom 2. The following "minimal" dimensions on the 
2-disk E = D occur: 

(i) dim V(D, 0) = 1 

(ii) dim V(D,1) = 

(iii) dim V(D, 2) = 1 

(iv) dim V(D, 3) = 1 

(v) dim V(D, 4) takes the minimal value consistent 
with (i)-(iv) (to be computed below). 

We make the further technical assumption that the 
constants a, b and c in Figure[S]are neither nor infinity. 

Axiom 2 (ii) is the "no tadpole" axiom which says that 
although figure ^b) is in Tt it has high energy. There is 
no low energy manifold V whatsoever when the boundary 
condition only allows tadpoles. If one thinks in terms of 

1 + f dimensional physics, 2 (ii) merely says the obvious: 
a single particle should not come out of the vacuum. If 
we nevertheless persisted in making dim (V(D, 1)) = 1 
we would admit the very boring case in which for all 
E, dim(V(E)) = 1 and V is spanned by the constant 
function on nets. 

Similarly 2(i) and 2 (iii) are required in a 1 + 1 dimen- 
sional (unitary tensor category) context. Also, if either 
dimension is then all V have dimension 0, by gluing 
formulae, so the entire theory collapses. 

Axiom 2(iv) does represent a choice. If we instead 
said dim (V(D,3)) = we would forbid our nets to 
branch. Here we know Z2 gauge theory (i.e. the toric 
code) and the doubled semion theory — both abelian 
- can arise. Possibly higher doubled SU(2) k Chern- 
Simons theories might also arise from loop models, but 
entropy arguments^ show that their ground state wave 
functions cannot be a simple Gibbs factor per \oopil. So 

2 (iii) is not inevitable but represents our decision to set 
up whatever microscopies are necessary to build (H, H) 
with branched nets occuring in V, i.e. to build a string 
net model. We now derive: 

Theorem 1 There is a unique theory V(E) compatible 
with 1 and 2 above. It satisfies dim (V(D, 4)) = 2 and is 
DFib, the doubled fibonacci category. 

Proof: By "Axiom 2" there are nonzero constants a, 6, c € 
C such that the conditions in Figure 00 hold. 

Now consider the 4-point space V(D,4). If its dimen- 
sion is < 4 there must be relations among the 4 nets in 
the expression below (fig. El a)), each thought of as an 
evaluation of the functional ^g.s. on V(D,4). We find 
conditions on the coefficients by joining various outputs 
of R. This amounts to calculating consequent relations 
in the 1-dimensional spaces V(D, 2) and V(D,3) which 
are implied by R. We work out, in part b) of fig. EH the 
implication of joining the upper outputs of R by an arc. 
By such arguments, we also have the three relations in 
fig. HlOfl . Eliminating x and y, x = —ch—bi, y = —bh — ci, 
we obtain 

(b - c- ab)h+ (-6 - ac)i = 

\-b - ac)i + (b - c - ab)h = 0. (5) 
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a (f> 




FIG. 8: We have stopped drawing the disk, but all diagrams 
above are nets in a 2-disk D with the endpoints on D. Also 
we abuse notation to allow a net to also represent the ground 
state evaluation, ^'g.s.(net). 



Possible relations which can be imposed on string nets 
through the Hamiltonian correspond to solutions of this 
linear system. There can be at most two linearly inde- 
pendent solutions, and this case occurs exactly when the 
coefficients vanish: b = c + ab, b = —ac so — ac = c+ (a 2 )c 
or a 2 = a + 1. We already see the golden ratio r emerge: 
a = t or -r -1 , t = 

In order to calculate the 6j-symbol, we set h = 0, i = 1 
in R to get the equation in fig. llll a). Unitarity requires 



1 so a — t and b 



yfr where <p is 

an irrelevant phase (associated with identifying the sim- 
plest 3-point diagram with some unit vector in V(D, 3)) 
which we set to 1. 



a) 



a) 



+ Y 



b) 



o 



bh 



+ + x 



+ ay 



bh + x + ay 



FIG. 9: If V(D, 4) < 4, then there must be some coefficients 
h,i,x,y such that the linear combination on the right-hand- 
side of (a) vanishes. By embedding these pictures within 
larger ones in such a way that the endpoints are connected 
as shown in (b) and using fig. |H| we find relations satisfied by 
h,i,x,y. 



=> ib + ax 



=> ch + ib 



=> bh + ci + y 



FIG. 10: Three other relations obtained by connecting the 
endpoints of R. 



+ b 



b) 



-1/2 



FIG. If: The 6j symbols can be obtained from the h — 0, i 
1 and h = l,i — a -1 relations. 



Similarly setting h = 1, i = a -1 (and hence x = 0) we 
find the second row of the 6j-symbol (i.e. the F matrix) 
in fig-E]b) and thus: 



F 



-1/2 



(6) 



Thus, fixing V (D, 4) to have minimal dimension, 2, we 
find our constants specified: a — T,b — r 1 / 2 , c = — t -1 / 2 , 
and the F-matrix as well. 

What we have obtained is the Turaev-Viro (or "dou- 
bled" ) version of the unitary Fibonacci fusion category. 
From these rules - a, b, c, and F - all nets G on a sphere 
can be evaluated to a scalar (G) which is the "golden" 
quantum invariant. Similarly, with this data the entire 
unitary modular functor is specified on all surfaces with 
or without boundary: V = DFib. 

Remark: If V{D,A) is allowed to have dimension 3, a 
generic solution, the Yamada polynomial, exists. It can 
be truncated to doubled SO(3)k - modular functors for k 
odd > 3 by imposing the correct dimension restriction on 
V(D, 2k). If V(D, 4) is allowed to be 4-dimensional, a re- 
lation in V(D, 5) realizes Kupperberg's G2-spider which 
presumably admits further specializing relations which 
generate Gi level k topological quantum field theories 
(TQFTs). 
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IV. DERIVING THE PROPERTIES OF DFIB 
FROM V(E,n) 

In the remaining three sections of this paper, we will 
discuss the properties of the ground states V(E, n) C 
H(T,,n), thereby obtaining physical properties of the 
topological phase DFib. DFib is the the product of two 
copies of opposite chirality of the Fibonacci theory, Fib, 
which the simplest 2 + 1 dimensional TQFT theory with 
nonabelian braiding rules. (Fib is also the simplest uni- 
versal theory^.) It arises as the "even sub-theory" of 
SU(2) S (i.e integer spins only) or 5(7(3)2 and also di- 
rectly from G(2)\. 

Fib has one non-trivial particle r with fusion rule: 



r = 1 



(7) 



Of course the quantum dimension of r is 1+ 2 = T - Q- n 
a slight abuse of notation, we use r to denote both the 
particle and its quantum dimension, the golden ratio.) 

A discrete manifestation of this quantum dimension is 
that Fib(S' 2 , n + 2), the Hilbert space for n + 2 r-particles 
at fixed position on the 2-sphere, is Fib(n), the nth Fi- 
bonacci number. (Proof: fuse two of the particles: the 
result will be either n + 1 or n r's on S 2 depending on the 
fusion process outcome. This yields the famous recursion 
formula: dim Fib(S' 2 , n + 2) = dim Fib(S" 2 , n + 1)+ dim 
Fib(S 2 ,n) which defines Fibonacci numbers.) This gives 
the exponential growth referred to in the introduction. 

Fib is a chiral theory with the following parameters: 



1 



V^+2 



1 r 

r -1 



= 37ri/10 



(8) 



(9) 



disks {Si}. According to Reshitikin-Turae\*i&, labelings 
(by 1 or t) of the bonds of G consistent with Fib fusion 
rules span the Hilbert space for a very high genus surface 
E^ which is the boundary of G, (G is a 3-D thicken- 
ing of the net G). It can be seen that the fixed space 
under the F-matrix action on G-labelings can be ob- 
tained from projectors associated to {Si}. These disks 
determine projectors, or plaquette operators, onto the 
trivial particle type along a collection of "longitudes" 
on E^ (where Si intersects E^). Another way to im- 
plement these plaquette operators is to add {Si} (thick- 
enings of {Si} to G). GlJ{<5i} is homeomorphic to a 
product E x 7, surface cross interval. Adding the pla- 
quet operators has cut Fib(E' t ) down to Fib(<9(E x 7)) = 
Fib(E ]J E) = Fib*(E) <g> Fib(E) = DFib(E). 

From this point of view the double arises from the fact 
that surface x 7 has both an "inner" and "outer" bound- 
ary. It is nontrivial to align the various structures (e.g. 
particles) of DFib from the two perspectives: one as a 
theory of trivalent graphs on a surface (Section III) and 
the other as a tensor product of a chiral theory and its 
dual. For this reason we are not content to merely state 
that the particle content of DFib is l®l,l®r, r®l, and 
t ®t. Rather, we will give a direct string representation 
for these particles shortly. 

Remark: The expression of DFib (and more genreally, 
Turaev-Viro theories) through commuting local projec- 
tors was known to Kitaev and Kuperberg and made ex- 
plicit for DFib in ref. 1. A conceptual understanding of 
the commutation relations is readily at hand from the 
preceding construction of E x I. The "longitudes" on 
which we apply plaquette projectors are disjoint and thus 
commuting. The fusion rules that are enforced at disjoint 
vertices also commute. Finally vertex and plaquet terms 
commute because rules are preserved under the addition 
of an additional ("passive") particle trajectory, labeled d 
in fig. El 



F 




T T 

r-V2 



-1/2 
1 



exp(4 pi i / 5) 



(10) 




FIG. 12: The addition of a passive arc to the (piece of a) 
string net on the right preserves the rules satisfied by nets. 




exp(2 pi i / 5) 



The theory V constructed in section 777 is isomor- 
phic to Fib* ® Fib = End(Fib). Very briefly we explain 
this connection in the context of a closed surface E (vi- 
sualize E = torus). Let G be a fine net on E which 
"fills it" in the sense that all complementary regions are 



We conclude this section by finding the 4 irreducible 
representations of the DFib annulus category. This is 
the linear G*-category whose objects are finite point sets 
(boundary data) on a circle and whose morphisms are 
formal combinations of nets in the annulus, which obey 
the linear rules a, b, c, and F, and which mediate between 
the boundary data. 

The four irreducible category (or "Algebroid") repre- 
sentations are detected as idempotents in Aq^ and A^i, 
the algebra under stacking of nets in annuli with either 



trivial or 1-point boundary data. A table which orga- 
nizes the results and compares back to the Fib* (g> Fib 
picture is given in the table below. The entries show the 
dimensions of the Hilbert space of (formal combinations 
of) nets on an annulus which start on the inner bound- 
ary with a given boundary condition (horizontal axis) 
and terminate near the outer boundary with a copy of a 
given idempotent (vertical axis). 



irreps (idempotents) 



boundary conditions 
12 



e x = 1 ® 1 

e 2 = t <g) t 
e 3 = 1 <g> t 
e 4 = 



1 



1 > 

1 1 > 

1 > 

1 > 



> 




FIG. 14: a) 1 G b) T G Ai,i, c) T~ 1 G A 1A 



(hint: apply F-matrix 
here) 



L = 




FIG. 15: A compact representation of L = r 1/2 l+r 3 < /2 (T+ 
T _1 ), as may be seen by applying the F-matrix where indi- 
cated. 



Let us start by finding the idempotents for A fi. A 0:0 
has the empty net as its identity and is generated by a 
single ring R. The a, b, c, F rules show: R 2 — 1 + R so 



b) 





FIG. 13: a) 1 G A , . b) R G A , . 



A 0;0 <* C[R]/(R 2 = 
More generally, 
C[i]/P(j),P(i) = 
tinct, are given by: 



1 + R). 

the idempotents 
: (x - ai) . .. (x - a k ) 



in the algebra 
all roots dis- 



(X - di) . . . (x - di) . . . (x - Ofc) 

- ai) . . . (a* - aj) . . . (ctj - Ofc) 



(11) 



We get 



ei 
e2 



1 + rii 
t + 2 



( = 1 (£> 1, the trivial particle) (12) 

( = T(g)T). (13) 



t = — t . Usinj 
bra A\ t \ is seen to 
and T 11 (where T 
an element L = t 
tors in a category 
the the elementary 
ful identities in Fig 
alently represented 



; the a,b,c, and F rules the alge- 
be generated by the identity 1, T, 
is defined in Figure I14|l . There is 
~ 1 / 2 1 + t- z ' 2 {T + T~ x ) which fee- 
sense through Aq.o- In fact, using 
rules in Figure OD we derive the use- 
ure 1161 and we find that L is equiv- 
as shown in Figure El It follows 



that to find a new representation of the annulus cat- 
egory we should look for the idempotents in A± t i/L. 
Ai.i/L = C[T]/{T 2 + tT + 1 = 0}. So, using again 
and a little manipulation we find idempotents 



e3 



and 



(-V?=3) 



t + ( T+V r^ ) i 



e 4 = 



(14) 



(15) 



in the quotient algebra. 





= - x 



FIG. 16: The relations depicted above can be obtained by 
applying the a, b, c, and F rules. 



Also L 2 = (t 1 ! 2 + t-3/2)L so e L = L/t 1 



/2 . 



.-3/2 

so again using the identities in Figure 1161 we can find 

e3 = e 3 (l - e L ) and e 4 = e 4 (l - e L ): 



e 3 = e 3 



2 + Vr^3 



2(7.1/2+7.-3/2)^ 



(16) 
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G4 = e<4 



(r + 2) - 



2( r i/2 + r -3/2) v ^— 3 



(17) 









p 2 









1/d 



U 

n 



CHROMATIC POLYNOMIAL AND 
YAMADA POLYNOMIAL 



FIG. 18: Jones- Wenzl projector. 



The chromatic polynomial Xg(^) °f a graph G at the 
positive integer k counts the number of fc-colorings of the 
vertices of the graph (so that no two vertices connected 
by a bond are given the same color). \ obeys the famous 
"delete-contract" recursion relation: 



(18) 



This relation can be depicted graphically as shown in fig- 
ure l|17|) . in which we have suppressed \ (as we have con- 
sistently suppressed the wave function and written the 
relation out in terms of its pictorial argument). This is a 



a) 



V 

A 



b) 



an idempotent familiar from the study of the Temperley- 
Lieb algebra TLj- Given G, ((G)) d is defined by labeling 
every arc of G by 2 as in l |18| l and then expanding to 
a weighted superposition of multi-loops. In each term, 
each loop contributes a numerical factor of d, 



((G)) d = 5>oeff.)d# lo °P s 



(20) 



For example, if G is a graph shaped like the Greek let- 
ter 6, we have the result shown in Figure El As this ex- 
ample shows, the Yamada polynomial is actually a poly- 
nomial in d, e? -1 . 



3/d 




+ 3/d 




-1/d 




FIG. 17: (a) Graphical depiction of the "delete-contract" re- 
cursion relation. This relation together with the two depicted 
graphically in (b) and (c) and multiplicativity under disjoint 
union completely determine X<g(fc). 



3 -1 -1 

= d - 3d + 3d - d 

3 -1 

= d - 3d + 2d 
FIG. 19: The Yamada polynomial for the 9 graph. 



local relation; there is no significance to the bit of G near 
the active edge e, which in the drawing is represented, 
purely for illustrative purposes, by three half-edges on 
top and bottom, x(^) is completely fixed by 1)18(1 and 
the following conditions: that x vanish on any graph in 
which the two ends of a single bond are joined to the 
same vertex, that Xsingie pt.(fc) = k, and multiplicativity 
under disjoint union: 



(19) 



We will be interested in Xq evaluated at noninteger val- 
ues as well as integral ones. 

We now turn to the Yamada polynomial defined for a 
net (trivalent graph) G lying in the plane (or 2-sphere); 
we denote it by ((G)) d , where d is the variable. To de- 
fine ((G)) d , recall the 2-strand Jones- Wenzl projector, 



Theorem 1 If G is a net in the 2-sphere and G is the 
dual graph then: 

((G)) d = d- v ^ XG (d 2 ) (21) 
V = number of vertices of G — number of faces of G. 

Proof: The above procedure for turning a net into a su- 
perposition of multi-loops may be generalized by declar- 
ing two local rules shown in fig. 1201 The first rule says 
a point is replaced by a circle with a possible numerical 
weight (vertex fugacity) and each arc is replaced by two 
lines with some general recoupling. We express the fact 
that this recoupling is still a variable to be solved for by 
Figure [21 Our goal is to find suitable values for p, u, 
and v so that Xg(^) comes out related to ((G)) d . The 
correct k will turn out to be d 2 . 



dual graph 



vertex 



multi-loops 




edge 



<0 



FIG. 20: Local rules for turning the dual graph of a net into 
a superposition of multi-loops. The ellipse represents an un- 
known combination 



CI) 



u 
n 



FIG. 21: A recoupling rule by which the dual of a trivalent 
graph can be converted into multi-loops. We will choose u and 
v (and also p) so that the procedure for making a trivalent 
graph into a multi-loop relates Xg(^) an d ({G)) d 




FIG. 22: Applying the rules in fig. I2UI to the graph on the 
left-hand side of the figure above yields the picture on the 
right. 



Let us look near a typical edge e of G and expand it 
using the rule in fig. [221 

With this expansion and fig. [21 we obtain fig. 123 
Translating this back into graphs yields fig. |2] Also, 




+ u p 




+ v p 



FIG. 23: Using the recoupling rule in Figure PHI we can sim- 
plify the picture on the right-hand-side of fig. 1221 

since graphs with an edge that connects a vertex to it- 
self evaluate to zero, we have fig. 1251 implying v = —ud. 
Also, k = pd. 

Comparing fig. [^] with the chromatic relation (|18fl we 
find: v = — l,u = 1/p. Because we have v = —ud, we 
obtain u = 1/d so d = p. Finally, since k = pd, k = d 2 . 
The mystery combination turns out to be a "sideways" 



up 



A 



FIG. 24: A relation that must be satisfied by u, v, and p. 



= 




ud f^3 




=> 



v 



ud 



FIG. 25: A graph with an edge that connects a vertex to itself 
evaluates to zero, from which we deduce a relation between u 
and v. 



P2 as shown in fig. [213 



CI) = -Vd 



^0 



FIG. 26: Solving for u,v,p, we see that the recoupling rule is 
just a "sideways" P2. 

To complete the proof, it remains to see the global 
geometry of how these P2 , s hook together. We claim 
that they lie along (doubled) dual graph edges. It suffices 
to examine an example. We take G and G to be the 
complete graph C4 shown in fig. [23 The factor of p v ^ 
appearing with ((G)) d has been put on the right-hand- 
side in the theorem statement. 

It is well known that the Yamada polynomial and the 
unitary invariant (G) T of Section III are closely related. 
Specifically, when d — t, we should modify the Ya- 
mada polynomial ((G)) T by a vertex fugacity to obtain 
(G) T . This can be seen by checking that both quanti- 
ties satisfy the elementary rules of Section III that de- 
termine them uniquely. From our example, shown in 
1191 we find ((d)) T = t _1 while for the unitary theory 
(0) T = ab = t 3 / 2 . To convert from the unnormalized 
(Kauffman) theory to the unitary theory (compare to 
ref. [l|) one must multiply in a factor of t 5 / 4 for each ver- 
tex. Thus: 



(G) T = (r^) v(G] r- V{6) Xa(r 2 ) 



(22) 



Using the Euler relation F(G) + V(G) - E(G) and the 
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net 



FIG. 27: Applying the rules to a complete graph, 
fact that G is a triangulation, implying E — |_F, we find: 
F(G) = 2V(G) - 4 (23) 

Therefore: 

(G) t = t-M^ X6 (t + 1) (24) 

We can now use Tutte's "golden ratio theorem" ([T] and 
[L]): For a planar triangulation G: 

( X6 (r + l)) 2 r 3V (r + 2)(r- 10 ) = Xd (r + 2) (25) 

This formula allowed Tutte to conclude that the r.h.s. is 
positive, creating a curious analogue (and precursor) to 
the 4-color theorem. (Neither result has been shown to 
imply the other.) 

Now square (|2~1|) and substitute into (|2*5|l ; the result is 
a remarkable formula: 

((G) t ) 2 = ^Xg(t + 2) (26) 

We have just proved the formula when G is a triangula- 
tion; in fact it holds more generally whenever G is a net. 
To establish this it suffices to check the formula when G 
is a single loop and to observe that the formula behaves 
well under disjoint union of disconnected components of 
G: the l.h.s. is obviously multiplicative and, it turns out, 
so is the r.h.s. The reason is that disjoint union of G\ 
and Gi corresponds to a 1-point union G\ \j G2 , and with 
the factor of 1/ (r+2), the chromatic polynomial becomes 
multiplicative under 1-point unions. So we have proved: 

Theorem 1 Let G be a net in the plane or 2- sphere (pos- 
sibly disconnected and possibly with circle components). 
Then 

(G) T 2 =(^)x 6 (- + 2) (27) 



Note: We would like to thank P. Fendley and E. Frad- 
kin for a discussion of this identity. In their papes^, a 
non-unitary normalization for (G) T led to a vertex fu- 
gacity on the right-hand-side of l|27|l . thereby obscuring 
the simplicity of (|27() and the direct connection to the 
Potts model. 



VI. STRING NET WAVEFUNCTIONS AND 
THE POTTS MODEL 

Let us review the high- and low-temperature expan- 
sions of the Potts model. We begin by assuming a lattice 
with Q "spin values" cr; = 0, 1, . . . , Q — 1 at each site i. 
The partition function Z(Q) is defined by: 

Z = ^>xp (-(3 \-jJ^S ai , a \\ (28) 

= EIT ex p(^^) ( 29 ) 

where the sum is over spin state configurations. Setting 
7 = e' 3,7 — 1 and expanding in powers of 7 we obtain the 
high-temperature expansion: 

z = Y J Q c i h (30) 

where the sum is over bond configurations; c is the num- 
ber of clusters, and b is the number of bonds. From 
now on, we consider the ferromagnetic case J = 1. Note 
that the last sum is over the 2 b distinct subsets, not the 
q(# of sites) distinct colorings because a Fubini resum- 
mation has taken place. Also, note that isolated sites 
count as clusters in II30JI . This is the Fortuin-Kateleyn 
representationii. It is known that for < Q < 4 the 
model is critical precisely at its self-dual point, 7 = \[Q- 
The high temperature expansion has been used2*£ to 
study loop gases with ground state wavefunction whose 
amplitude is d L for some real number d 1 where L is the 
number of loops. The square of such a wave function 
can be interpreted as a Gibbs weight (d 2 ) L providing a 
"plasma analogy" between topological ground states and 
the statistical physics of loop gases. We can easily see 
that a loop gas with Gibbs weight d 2 per loop is critical 
if d < V2: 

(d 2 ) L = (d 2 ) c+c " = (d 2 ) 2c+b = (d 4 ) c (d 2 ) b (31) 

Here, c* is the number of dual clusters, i.e. the minimum 
number of occupied bonds which have to be cut in order 
to make each cluster tree-like (essentially the number of 
"voids" which are completely contained within clusters). 
The first equality follows from each loop being either the 
outer boundary of a cluster or a dual cluster. The second 
is due to the Euler relation c* = c + b + const. To map 
this squared wavefunction to the Potts model, we need 
d 4 = Q < 4; note the edge fugacity for this loop model 



11 



is automatically y/Q, placing the model at its self dual 
point. 

The low-temperature expansion of the Potts model is: 
Z(Q)=J2xdQ)(^ J ) L (32) 

G 

where the sum is over graphs G and L is the total length 
of the graph. We specialize to the case in which the dual 
lattice is trivalent (e.g. the triangular lattice, whose dual 
lattice is the honeycomb), so that the graphs will are all 
trivalent. At the critical point (Q < 4), this becomes 

using the condition for criticality e^ cJ — 1 = y/Q. 

The (unnormalized) isotopy invariant wave function W 
constructed in Section III satisfies \V(G) = (G) T . Hence 
the corresponding statistical physics (of equal time corre- 
lators) is that of the normalized probability distribution: 

Pr ° b(G) = Z(Q = r + 2,^ = 0) ^ ™ 

Note that (using 7$) this is the high temperature limit 
of the low temperature expansion. 

This situation reminds us of the Toric codei^, where 
the weight on loops may be obtained by setting Q = 2, 
(3 = (for Q = 2, branched nets have zero weight), 
so the toric code ground state wavefunction can also be 
understood in terms of the high-temperature limit of the 
low-temperature expansion of the Potts model. In both 
cases, we conjecture that for (3 < (3 C the wavefunction 
with isotopy invariance modified by: 

*(G) = (G) r :r lcnsth(G) (35) 

for 1 > x > l/(y/Q + 1) (w 0.345 when Q = r + 2, 
and » 0.466 when Q = 2) will be the ground state for 
some gapped Hamiltonian in the corresponding topolog- 
ical phase DFib (or Toric code). 

Numerical workiS already supports this conjecture in 
the Toric code case. Also, note that at [3 > (3 C (so x < 
l/(v / Q + 1)) an effective string tension prevents the nets 
from fluctuating. Thus, the system leaves the topological 
phase and enters a "geodesic phase" in which small nets 
dominate. Recoupling is now unlikely, so the necessary 
topological relations are not well enforced and geodesic 
continuation allows states on the torus to be fairly well 
guessed by measurement within a subdisk violating the 
disk axiom. Finally, (|27|) identifies the ground state of 
the gapped Levin- Weni Hamiltonian: \V(G) = (G) T as 
the (3 = end point of the conjectured family. 

In ref. @, it is conjectured that the unnormalized prob- 
ability distribution on nets (((G)) d ) 2 should somehow 
yield the same statistical physics as the Potts model at 
an "effective" Q c s satisfying: 

(Qcff-l) = (d 2 -l) 2 . (36) 



If we set d = t and adjust the G-vertex fugacity so as to 
replace ((G)) T with (G) T , as we did in Section V, then 
the conjecture is precisely verified atQ = r + 2,/3 = 0: 

(t + 2-1) = (t 2 -1) 2 (37) 

We close this section with an observation whose impor- 
tance, if any, we do not yet understand. The high temper- 
ature expansion weights loops while the low temperature 
expansion weights nets. But loops are nets, so we might 
ask: at what value of Q do the high and low temperature 
expansions weigh loops equally? In the low temperature 
expansion, loops have weights Q — l and in the high tem- 
perature expansion they have weight y/Q (from UTT)) . So 
loops are equally weighted when Q — 1 = y/Q, i.e. when 
Q = r 2 . 



VII. CONCLUSIONS 

DFib is among the simplest conceivable achiral particle 
theories. In some sense it rivals the toric code in simplic- 
ity (both have 4 particles), but is vastly richer (in fact, 
universal^) in its braiding. We have explored the path 
to this phase, taking a Hilbert space H of nets and an 
isotopy invariant Hamiltonian (possibly with a bond fu- 
gacity, see Section VI) H : H — > H as our starting point. 
The exploration has been combinatorial in Section III, 
algebraic in Section IV, and statistical in Section VI. 

In Section III, we saw that DFib emerges from mini- 
mizing certain degeneracies. This encourages us to be- 
lieve the phase will ultimately be found in nature - as 
nature abhors a degeneracy. 

In Section VI, we establish that the nets G with 
(squared) topological weighting ((G) T ) 2 (as usual, squar- 
ing the wavefunction to obtain a probability) are in a 
high temperature phase of the (r + 2)-state Potts model 
(above criticality). This is also encouraging: the classi- 
cal critical point looks as if it is the "plasma analogy" 
of a quantum critical point sitting at the entrance to the 
DFib phase. A parallel is explored between this situa- 
tion and the Q = 2 Potts critical point which serves as 
an entrance to the toric code phase. 

Unresolved is what, more precisely, is required of H : 
Ti — > Ti to be in the DFib phase. Is enforcement of the 
net G structure (encoded in the definition of Ti) plus 
strong dynamic fluctuation of G adequate? We do not 
know. 
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APPENDIX A: THE CHROMATIC 
POLYNOMIAL AND HARD HEXAGONS 

In this appendix we extend a theorem of Tutted which 
bounds the decay of the chromatic polynomial at t + 2 
of planar graphs. Specifically, we obtain a better sharp 
bound for graphs that consist of a large regions of hexag- 
onal lattice (see Figure 1281) . The ground state of the 




FIG. 28: large chunk of hexagonal lattice 



doubled Fibonacci theory DFib consists of a certain su- 
perposition of string net configurations on a hexagonal 
lattice. As in ref. we may construct from it a classi- 
cal statistical mechanical model of nets, with the Boltz- 
mann weight of each net equal to the norm squared of 
its ground state amplitude. As discussed in Section V, 
this statistical mechanical model is the (solvable) infinite- 
temperature limit of the Q = r + 2 Potts model. We are 
interested in the rough quantitative behavior of the am- 
plitudes for different types of graphs. In particular, we 
can consider a large chunk of hexagonal lattice, which is 
just the string net consisting of a large finite region whose 
bulk includes every available bond (see Figure!^. We 
can ask about how its amplitude scales with the size of 
the region. Now, the Boltzmann weight of a string net is 
(up to overall scaling) just the chromatic polynomial at 
t + 2 of the graph dual to the net, which can be inter- 
preted as the zero=temperature limit of the r + 2 anti- 
ferromagnetic Potts model defined on the graph dual to 
the net. So for the large chunk configuration, we are just 
interested in the behavior of the zero-temperature free 
energy of the r + 2 antifcrromagnetic Potts model on the 
hexagonal lattice. 

This model is criticaUi and is expected to be described 
by a conformal field theory. On general grounds^ we 
expect the free energy to scale like 

F = cqA + ciL + c 2 log L (Al) 

where c-i is universal and related to the central charge 
of the CFT and Co and c\ are non-universal (here A is 
the area, i.e. number of hexagons in the region, and L is 
the length of the boundary) . It turns out that the model 



is exactly solvable and we can obtain an exact analytic 
expression for cq. 

To do so, we first use the so-called shadow method^ to 
evaluate G T of a net. This method works as follows. We 
take the net to be located on the sphere and assume that 
the regions (or faces) that it bounds are simply connected 
and do not border themselves. Then the shadow method 
(applied to DFib) gives G T as a sum over black and white 
colorings of the faces 

G T = FcEc^Vc. (A2) 

colorings C 

The colors black and white are identified with the two 
particle types of the theory, and 

Fc= II d c(F) (A3) 

faces F 

e c = n q ^ e ^ ( m ) 

edges E 

V c = Yi Tetrahedron(C, V). ( A 5) 

vertices V 

Here dc(F) stands for the quantum dimension associated 
with the color of face F (1 for black, r for white). 0(C, E) 
is the theta graph along whose three edges run particle 
types associated with the edge E (always the nontrivial 
particle) and the two faces which E borders. Tetrahedron 
(C, V) is the tetrahedral graph with its six edges labeled 
by the particle types corresponding to the faces and edges 
adjoining V. 

Let v — # vertices, / = # faces, e = # edges. When 
applied to the net which consists of a large chunk of 
hexagonal lattice, the shadow method reduces to the 
following: the region containing the point at infinity is 
black, and the sum over colorings becomes a sum over col- 
orings in which there are no adjacent black hexagons (and 
no black hexagons adjacent to the black outside region). 
Let us first compute the weight of the all white coloring. 
We ignore boundary effects. We have first of all v = 2/ 
and e = 3/. Each edge contributes _1 = t~ 3 / 2 , each 
vertex Tetrahedron = — t, and each face the quantum 
dimension r. The total weight of the all white coloring 
is thus 

(t-*' 2 )* {-T) v {r) f = t-V' 2 . (A6) 

Now suppose we have a configuration with some black 
hexagons. Its weight is just that of the all white configu- 
ration multiplied by appropriate ratios of theta symbols, 
tetrahedron symbols, and quantum dimensions. More 
specifically, for each white hexagon that one turns into a 
black hexagon, one must multiply the edge contribution 
by (t/t 3 / 2 )~ 6 , the vertex contribution by (t 3 / 2 /(— t)) 6 , 
and the face contribution by t _1 . The product of these 
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is t 5 , so that the weight of each such coloring is just 
T 5#( black hexagons) timeg the weight of the a ll-white col- 
oring. The sum over colorings now just yields the critical 
Hard Hexagon model, whose free energy per vertex was 
obtained by Baxter— (eqn. 10). We thus obtain 



G 



f 



where 



'27(25 + 11^5)' 
250 



1/2 



From we surmise 

X (r + l) = G T r^ 5 / 4 " = (r- 3 Kc ) / . 



(A7) 



(A8) 



(A9) 



We note that this is sharper than Tutte's bound^ of 
const. r~f on x(r + 1): 



0.546 / < 0.618 / . 



(A10) 



Tutte's bound is of course more general in that it applies 
to any net configuration on the sphere (whose regions 
are simply connected and don't border on themselves). 
For completeness, we note from Ij26|) and IjAip that Co = 

T- 3 /V. 
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